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ABSTRACT 

If the dark energy equation of state is anisotropic, the expansion rate of the universe becomes direction- 
dependent at late times. We show that such models are not only cosmologically viable but that they 
could explain some of the observed anomalies in the CMB. The possible anisotropy can then be 
constrained by studying its effects on the luminosity distance-redshift relation inferred from several 
observations. A vector field action for dark energy is also presented as an example of such possibility. 
Subject headings: (cosmology:) cosmic microwave background; cosmology: miscellaneous; cosmology: 
observations; cosmology: theory; cosmology: large-scale structure of universe 
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1. INTRODUCTION 

It is not clear whether the large-angle anomalies in the 
observed cosmic microwave background a re of a cosmo- 
logical origin and not due to systematics jEriksen et al.l 
120041 lLand and MagueikJ2005l : ICopi et al.ll2006f ) Never- 
theless, it is tempting to associate the apparent statistical 
anisotropy with dark energy, since the anomalies occur 
at the largest scales, and these enter inside the horizon at 
the same epoch that the dark energy dominance begins. 

The paramount characteristic of dark energy is its neg- 
ative pressure. We examine then the possibility that this 
pressure varies with the direction. Then also the uni- 
versal acceleration becomes anisotropic, and one would 
indeed see otherwise unexpected effects at the smallest 
multipoles of CMB. At the background level we will find 
a similar CMB pattern as in an uni yerse which is el- 
lipsoi dal at the era of last scattering (|Campanelli et ail 
2006). We then show that the supernovae could be used 
to distinguish these different scenarios and to constrain 
the possible anisotropic properties of dark energy. 

There are several motivations for anisotropic models 
of dark energy. For instance, many dark energy mod- 
els are in principle compatible with the FLRW metric 
but exhibit anisotropic stresses at the perturbative level, 
including nonminimally couple d fields, viscous fluids 
and modified gravi t y models (IKoivisto and Motal 120061; 
Schimd et all 120061: UMota et aLM2007t iKoivisto & Motal 



2007allbt ). It is possible that a more accurate description 



of such models should take into account the anisotropic 
effects on the background expansion, which then breaks 
the statistical isotropy of perturbat ions too (as appears 
to have happened in the CMB (lEriksen et alt 12004 
lLand and Magueiiol [20051 : iCopi et al.ll2006D ). Also, con- 
sidering dark energy as an effective description of a 
backreaction of sizable inhomogeneities in dark m atter 
(j Buchertj|2007t I Behrend. Brown fc Robbers! l200l . the 
validity of a perfect fluid description could seem dubious. 

To have a general description of an anisotropic dark 
energy component, we consider a phenomenological pa- 
rameterization of dark energy in terms of its equa- 
tion of state (w) and two skewness parameters (5, 7) 



and include also a coupling term (Q) between dark en- 
ergy and a perfect fluid (dark matter). Previous stud- 
ies anisotropic dark energy have mainly considered the 
anisot ropic properties of the inhomogeneous perturba- 
tions ()Koivisto and Motal 120061 iBattve and Mossll2006f) 
whereas our approach here is to focus on a smooth cos- 
mology with the anisotropic pressure field. Whereas in a 
previous approach the anis otropy was only weakly con- 
strained (jMota et al.ll2007f ), we find that in the present 
description only a narrow parameter range can survive 
the observational tests. Similar anisotr opic inflation has 
been considered in the early universe (|Burd and Lidsevl 
fl99il iFordl |l989) . Therefore we focus on the present 
paper more on the dark energy era, where the qualita- 
tive differences are that matter cannot be neglected and 
that the cosmologies do not isotropize. As a proof of a 
concept, we write also an explicit field theory example 
where a vector field drives the anisotropic acceleration of 
the universe. 

An anisotropic expansion is not compatible with the 
Robertson- Walker (RW) metric. Hence we use the B(I) 
(Bianchi type I) metric which generalizes the flat RW 
metric and may be employed to obt ain limits on c osmo- 
logical skew pressures from the CMB (|Barrowl ll997) . The 
line-element of a B(I) universe is: 



ds 2 = -dt 2 



a 2 {t)dx 2 +b 2 (t)dy 2 +c 2 (t)dz 2 . (1) 



There are thus three scale factors, and consequently three 
expansion rates. In principle all these could be different, 
and in the limiting case that all of them are equal one 
recovers the RW case. It is useful to express the mean 
expansion rate as an average Hubble rate H 



(2) 



(where an overdot means derivative wrt t) and then 
the differences of the exp ansion rates a s the Hubble- 
normalized shear R and S (|Barrowlll997f ) 
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We consider a universe filled with a per- 
fect fluid having the energy-momentum tensor 
T {m)u = diag(-l,w m ,w m ,w m )p m , and dark energy, 
which we allow to have the most general energy- 
momentum tensor compatible with the metric (JTJ) 

T>* v = diag{-l,w,w + 3S,w + 3j)p. (4) 

The generalized Friedmann equation may be written as 

InG p m + p 



H 2 



3 1 - i (R 2 + S 2 - RS) ' 



(5) 



We let the two components also interact. The continuity 
equations are then 



and 



p + 3 



3H (1 



(l + w)H- 



QHp, 



(G) 



c 



QHp, (7) 



where Q determines the coupling. If it vanishes, it fol- 
lows that p m ~ (abc)- 1 -™™ and p ~ {abc)- 1 - w b 3S c 3 ~< . 
Defining x = i log (— g), where the metric determinant 
g = —abc, one notes that H — x. We will use x as our 
time variable rather than t. Derivative wrt x is denoted 
by star. We also define the dimensionless density frac- 
tions 

8irG p m _ p 



3 H 2 ' 



U : 



Pm + P 

Using R, S, and U as our dynamical variables, the system 
can finally be written as 



U* = U (U - 1) [7 (3 + R - 2S) + 

S(3~2R + S) + 3(w- w m )] ~ UQ, 



S* = - (9 - R 2 + RS - S 2 ) 
6 

x {S [U (5 + 7 + w — w rn ) + 

R* = I (9 - R 2 + RS - S 2 ) 
6 



(8) 



1] - 67C/} , 



x {R [U (S + 7 + w - w m ) + w m - 1] - 65U} . 

Note the coupling term Q appears only in the evolution 
equation for U. Never theless its presence can ch ange the 
dynamics completely (|Koivisto and Motall2008f ). 

2. VECTOR FIELD 

As a proof of concept, we present an explicit field the- 
ory for the anisotropically stressed dark energy. Consider 
the vector field action 



S 



d 4 x. 



-9 



R 



16ttG 



-iV^ - V + qL, 



(9) 



where the kinetic term involves -F M „ = d^A v — d v A^. The 
potential V and the possible couplings q are understood 
as functions of the quadratic term A 2 = A^A^ . The 
field equations are = 8ttG (T™ + T*) , where the 
energy-momentum tensor of the vector field follows by 
varying wrt the metric 



T*=F lta F2 + 2V'A ll A v - 



\ F <# F 



a/3 



+ V g» v . (10) 



Prime denotes derivative wrt A 2 . With the metric ((T|), 
the equation of motion for the Aq component of the field 



dictates that [q'L m + V']Aq — 0. Thus we confine to 
purely spatial vector fields. The spatial components of 
the energy-momentum tensor are, 



T 



iA 



■Adj + 2V'A i A j 




The off-diagonal terms should vanish. Thus one is re- 
stricted to consider only vectors which are parallel to one 
of the coordinate axis. There, however, could be many 
such fields. 

In the RW case all the diagonal components of 
T\j should be equal, which indeed requires three vec- 
tor fields, one in each coordinate direction, of ex- 
actly equal magnitudes, as t he triads of Armendariz- 
Picon (|ArmendarizT > icon 2004). With the metric (JTJ) , the 
only constraint is that every vector field should be along 
a coordinate axis. Any system of these vector fields is 
thus required to have an anisotropic equation of state, 
unless it reduces to the very special triad case. Let us 
therefore briefly look at the case of a single vector field 
with A\ = A3 = 0. Notice that, in analogy to system 
((8]), the coupling Q is now related to the function q. If 
we define the additional dimensionless parameters 



Y EE 



V 



Y,=2- 



H 2 ' 1 b 2 H 2 ' 

one can relate the present case to the general notation 
by noting that H 2 p = X 2 /2 + Y, and that if S = 0, then 

X 2 -2Y „ -X 2 - Yi 



w 



X 2 + 2Y' 



3.X 2 



3Y 



(11) 



We see that the anisotropy is naturally small if the field is 
either subdominant or near its minimum. Such models, 
even though perturbatively close to standard cosmology, 
would be excluded by imposing the RW symmetry. Vi- 
able models with large anisotropy do also exist as shown 
below. An extensive inves tigation of specific model s will 
be undertaken elsewhere (|Koivisto and M ota 2008). 

3. SCALING SOLUTIONS (IN THE AXISYMMETRIC CASE) 

To begin, we will assume for simplicity that 1) the 
perfect fluid is minimally coupled dark matter, w m = 
Q = 0, and further that 2) the skewness together with 
the equation of state of dark energy is constant, w = 8 = 
7 = (see (jKoivisto and M ota 2008) for an extensive 
analysis with these assumptions relaxed). In our model 
the universe is then initially (close to) isotropic, but is 
driven to anisotropic expansion by the skewness of dark 
energy. Under the additional simplifying assumption of 
axial symmetry (S = 7 = 0), in FIG. [1] we indicate the 
final possible stages of such universe, corresponding to 
the fixed points of the system ([5]). The universe will 
then end up in three possible scenarios: (I) The isotropic 
RW case corresponding to 

R = 0, U = 0. (12) 
(II) An anisotropic dark energy dominated solution with 
66 



R 



w + S - 1 



U = 1. 



Or (III) a scaling solution 
3 (w + 5) 



R 



25 



U 



-jw + S) 
38 2 — 2Sw — w 2 



(13) 



(14) 



3 
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U=0 (FRLW) 



< U < 1 



U = 1 



-1 



Fig. 1. — Asymptotic state from an Einstein-deSitter stage. The 
future fate depends on the dark energy properties w and <5 and is 
classified into three possibilities. If <5 + w > 0, the isotropically ex- 
panding dust domination (solution (1 1 -'li for which U = 0) continues 
forever. Otherwise, the universe will end up expanding anisotropi- 
cally and either dominated by dark energy (solution 1131 for which 
U = 1) or exhibiting a scaling property (solution IU4H for which 
<£/<!). 



Given that w is negative, this solution accelerates w e f / 



2 H 

3 H 

that in the R 



1 < ~, if w < 8 < 



<S< 



Note 



case scaling solutions could only be 
found for coupled components. 

Notice that within the RW universe it has been proven 
difficult to address the coincidence problem by finding a 
model entering from a matter-dominated scaling solution 
to an accelerating scaling solution. Allowing for the pres- 
ence of three expansion rates opens up the possibility of 
describing a universe entering from a perfect fluid dom- 
inated scaling to an anisotropically accelerating scaling 
era. This might eventually help to understand the coinci- 
dence problem, since then matter and dark energy would 
have had similar energy densities both in the past and in 
the future (in the past the dark en ergy fraction, if con- 
stant should not exceed about 1/lO lDoran and Robbers! 
(|200l V 

We now study the observational implications of mod- 
els with nonzero skewness, relaxing the assumption 2). 
The relevant case is then a universe entering from RW 
solution (fT2"|) to an anisotropically accelerating universe 
(which generalizes (fl"3|) if 7 ^ 0). 

4. CMB ANISOTROPY 

The B(I) background predicts a quadrupole pat- 
tern in the CMB. It is not clear how the fluctua- 
tions in the photon temperature distribution evolve 
when the background expansion becomes anisotropic, 
since the anisotropic sources, even if perturbatively 
small, couple to the per turbations at the first order 
(|Koivisto and M ota 2008|). It seems inevitable that sta- 
tistically anisotropic features would then be, in prin- 
ciple, present in the whole spectrum of fluctuations 
(though in our case such features would be confined 
to larger scales as the anisotropy of the expansion be- 



comes important only recently). Sever al large-angle 
anomalies have been reported in the CMB (Eriksen et al.l 
12004 lLand and Magueiidl2005l ; ICopi et al.ll2006h . It has 
been demonstrated t hat CMB spectrum of the tilted 
B(VIIh) background rtHawkmdll968tlBarrow et alJll985t 
iPontzen and Challinorl l2007fl would hav e potential to 
model these anomalies (paffe et al.l l2005h were it not in 
conflict with a dark e nergy dominated universe at the 
smaller angular scales paffe et al.ll2006D . 

The B(I) model reproduces the predictions of the 
concordance model at small scales, while featuring 
anomalies at large angles, since the (originally statisti- 
cally isotropic) CMB field experiences an (statistically) 
anisotropic integrated Sachs- Wolfe effect in the nowa- 
days ellipsoidal universe. Furthermore, it is plausible 
that these anomalies can match the observed ones, since 
a B (Villi) model can be effectively described as a B(I) 
model with an additional anisotropic energy source 1 

As a first step, we check whether the background 
|T]) is compatible with the available cosmological 
data. There are a number of calculations of CMB 
anisotropies in gener al Bianchi uniyerses already avail- 
able in the literature (lHawking]|1968t iBarrow et al.lll985t 
IPontzen and Challinorll2007f) . Here we consider the dark 
energy created quadropole and the effects to the lumi- 
nosity distance. 

By considering the geodesic equation for photons, one 
can derive an equation for the redshift z of a photon 
arriving from the a direction p 



l + z(p) = -Jl 
a V 

in terms of the eccentricities 



= (-) -1. 

- C/ 



(15) 



(16) 



Note that the scale factors and eccentrities here are eval- 
uated at the time of last scattering in the case that the 
scale factors are all normalized to unity today. If T* is 
the temperature at decoupling, the temperature field is 
given by T(p) = T*/(l + z(f>)), and its spatial average is 
AttT = J d£lpT(p). The anisotropy field is then 



f 



1 



T(P) 
T 



(17) 



The coefficients in the spherical expansion of this 
anisotropy field are called ag rn , and due to orthogonality 
of spherical harmonics Yt m are given by 



= / dn, 



ST(p) 

> 

T 



Y, 



The multipole spectrum is described by 



1 £(£ + 1) 
\| 2tt (2£+l) 



E 



\O-lr, 



(18) 



(19) 



Expanding the redshifts (TT5)l in powers of the eccen- 
tricities (fl~6"|) one notes that the ai m will be real (since 
there is only even dependence on the polar angle in the 

1 One may obtain the B(VIIO) model from B(I) by adding 
anisotropic curvature. To obtain B (Villi) model one should add 
also an isotropic curvature similar to such present in B(V) models. 
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anisotropy field and the imaginary parts of the e ml ^ in- 
tegrate to zero), and that for all odd £ the ai m will 
vanish (since only even powers of the azimuthal an- 
gle appear in the expansion). To first order in e. 



T = aT* [1 — g (ef + e^)] , and in addition to the monopole, 
there is only the quadrupole 



O20 : 



2ei 



a 2 i = 02-1 = 0, 



(20) 



0,22 — 0,2-2 — — 



30 ^ 



Qi = 



5V3 



The observed value o f Q2 is lower than t he concor- 
dance model predicts (jHinshaw et al.1 [2007). It has 
been suggested in previous works that this discrep- 
ancy cou ld be explained by an ellipsoidality of the 
unive rse (|Campanelli et ail 120061 : iJimenez and Marotol 
2007). This would require that the anisotropy of the 
background is suitably oriented with respect to the in- 
trinsic quadrupole and cancels its power to a sufficient 
amount. Too large anisotropy would of course only make 
the situation worse regardless of the orientation. De- 
pending on the cosmological model, one should have 
Qi i$ 2 ■ 10~ 5 to be consistent with observations tak- 
ing into account the cosmic variance. The constraints 
this implies on the skewness of dark energy are very 
tight. However, we remark that in more general mod- 
els, in particular with time-varying S and 7, one could 
allow more anisotropy. It is in principle possible for ar- 
bitrarily anisotropic expansion to escape detection from 
CMB (considering only the effects from the background) , 
as long as the expansion rates evolve in such a way 
that e z — e y = 0. In other words, the (background) 
quadrupole vanishes, if each scale factor has expanded - 
no matter how anisotropically - the same amount since 
the last scattering. An example of such a scenario, de- 
rived from the action ((9]), is shown in FIG. [2] 

5. SNIA LUMINOSITIES 

The luminosity-redshift relationship of the SNIa data 
could be used to probe the possible anisotropies in the 
expansion history. This is a useful complementary probe 
since these objects are observed at the z < 2 region, 
whereas CMB comes from much further away at z ~ 
1100. The luminosity dista nce at the redshift z in the 
direction p is now given by (|Koivisto and M ota 2008) 



d L (z,p) = (1 + z) 



dt 



p 2 a 2 + plb 2 



p2 c 2 



(21) 



To test this prediction with the data, we apply the for- 
mula (|15p for each observed redshift of a supernova and 
match its luminosity distance inferred from the obser- 
vation to the one computed from (|2ip . In addition, we 
have to take into account also the angular coordinates 
of each individual supernovae in the sky to fix p for 
each object. In o ur analysis we use the GOLD data set 
(R iess et al1l2006l) . which consists of five subsets of data 2 . 
We marginalize over the directions in the sky and over 
the present value of the Hubble constant. The results are 
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Fig. 2. — A minimally coupled vector field satisfying the 
quadrupole constraint. The solid (black) line is the dark energy 
density fraction U, the dashed (red) line is the effective equation 
of state of the universe w e ff, and the dash-dotted (blue) line de- 
scribes the evolution of eccentricity, E = 5Q0ey . The potential is a 
double-power law V = mA 2 + Ayl -4 . The dynamics of the field is 
such that though there are significant anisotropies, the eccentricity 
at the present is close to zero. This may be achieved with different 
power-law potentials, but requires fine-tuning of the mass scales m 
and A. 



summarized in FIGs. [3]and|4] The best-fit anisotropic 
models are only slightly preferred over the ACDM, the 
difference being A^ 2 ~ 1. Because of the additional pa- 
rameters in the anistropic models, the reasonable inter- 
pretation of these statistics is that the SNIa data favors 
isotropic expansion. 
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Fig. 3. — Limits on the (constant) skewness parameters of dark 
energy arising from the SNIa data, when when f2 m = 0.3 and 
w = —1. The contours correspond to 68.3, 90, 95.4 and 99 percent 
confidence limits. 



SNIa Angular coo rdinates of each of the 182 GOLD s uper- 
novac can be found in URiess et al. 2006; Asticr et al. 200D and 



|http://cfa- www. harvar d.edu/ps/lists/Supernovae. html 



The SNIa data constrains the skewness parameters 
much looser than the CMB quadrupole, if they are con- 
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Fig. 4. — Constraints arising from the SNIa in the axisymmetric 
case 7 = 0. Inside the darker isosurface, the fit is as good as in 
the ACDM model, \ 2 < 158. Inside the lighter isosurface, one 
has A\ 2 < 8.02. One notes that larger skewness 8 would typically 
be compatible with the SNIa data for phantom equations of state 
w < — 1 and large matter densities. 



a way that it appears to us anomalous at the largest 
scales 3 . The future SNIa data, with considerably im- 
proved error bars ( e.g. from the SNAP experiment, 
|http://snap.lbl.gov/| ) on cIl(z,p) might be used to rule 
out this possibility, and to distinguish whether the pos- 
sible statistical anisotropy was already there at last scat- 
tering or whether it is due to dark energy. 



6. CONCLUSIONS 

In conclusion, dark energy with an anisotropic equa- 
tion of state might be the culprit for both the cosmic 
acceleration and the large- angle anomalies in the CMB. 
This might also be the key to understand the coincidence 
problem. The present SNIa data allows anisotropic ac- 
celeration, but SNAP could set things straight about the 
skewness of dark energy and so of its nature. Such pos- 
sibility would open a completely new window not only 
on the nature of the CMB anomalies but also into high 
energy physics models beyond the usual isotropic candi- 
dates of dark energy such as scalar fields or the cosmo- 
logical constant. 



stant. However, in specific models with time-evolving S 
and 7, the constraints from CMB and from SNIa can 
be of comparable magnitude and allow anisotropy to 
an interesting degree. This means that the even if the 
CMB formed isotropically at early time, it could be dis- 
torted by the acceleration of the later universe in such 

3 Generating these effects at low redshift has an advantage that 
it relaxes co nstraints which would otherw ise come from the CMB 
polarization (Pontzcn and Challinor 2007) and could be strong for 
a given temperature anisotropy in isotropizing models because of 
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